Low- Metallicity Star Formation: From the First Stars to Dwarf 
Galaxies 

Proceedings lAU Symposium No. 255, 2008 © 2008 International Astronomical Union 

L.K. Hunt, S. Madden & R. Schneider, eds. DOI: 00.0000/XOOOOOOOOOOOOOOOX 

The Star Formation History of the GRB 
050730 Host Galaxy 

Francesco Calura^ 

^ Dipartimento di Astronomia, Universita' di Trieste, via G. B. Tiepolo 11, 34131 TRIESTE - 

ITALY 

email: fcalura@oats.inaf.it 

Abstract. The long GRB 050730 observed at redshift z ~ 4 allowed the determination of the 
elemental abundances for a set of different chemical elements. We use detailed chemical evolution 
models taking into account also dust production to constrain the star formation history of the 
host galaxy of this long GRB. For the host galaxy of GRB 050730, we derive also some dust- 
related quantities and the the specific star formation rate, namely the star formation rate per 
unit stellar mass. We copare the properties of the GRB host galaxy with the ones of Quasar 
Damped Lyman Alpha absorbers. 
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1. Introduction 

Gamma-ray bursts (GRBs) afterglows provide an insight into the intersteUar medium 
(ISM) of galaxies during the earliest stages of their evolution. In several cases, thanks to 
the GRB afterglows it has been possible to study the dust content, the star formation 
rates and the stellar mass of the GRB host (Bloom et al. 1998; Savaglio et al. 2008). 
In a few cases, the determination of their chemical abundance pattern has been possible 
(Savaglio et al. 2003; Vreeswijk et al. (2004); Prochaska et al. 2007). 
Chen et al. (2005) reported on the chemical abundances for the damped Lyman Alpha 
system (DLA) associated with the host galaxy of GRB 050730. Their analysis showed 
that this gas was metal poor with modest depletion. These results were subsequently 
expanded and tabulated by Prochaska et al. (2007) (hereafter P07). 
In this paper, we aim at determining the star formation history of a GRB host galaxy by 
studying the chemical abundances measured in the afterglow spectrum of GRB 050730. 
For this purpose, we use a detailed chemical evolution model. Our aim is to constrain 
the star formation rate and the age of the host galaxy of GRB 050730, and possibly to 
expand our study to other systems. 

The plan of this paper is as follows. In Section 2, we briefly introduce the chemical 
evolution model. In Section 3, we present and discuss our main results. 

2. Chemical evolution models including dust 

The chemical evolution model used here to derive the star formation history of the 
GRB 050730 host galaxy is similar to that developed by Bradamante et al. (1998) for 
dwarf irregular galaxies. Since some chemical species studied in this paper are refractory 
(Fe, Ni), in the chemical evolution model we include also a treatment of dust production 
and destruction, based on the work by Calura, Pipino & Matteucci (2008, hereinafter 
CPM08). 
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We use the model for dwarf irregular galaxies since several observational investigations 
provided strong evidence that most of the GRBs originate in gas rich, star-forming sub- 
luminous {L < L^) galaxies with relatively low metallicities {Z < Zq) (Bloom et al. 
1998; Prochaska et al. 2004). 

The dwarf galaxy is assumed to form by means of a continuous infall of pristine gas until 
a mass Mtot is accumulated. The evolution of dwarf irregular galaxies is characterized 
by a continuous star formation history, characterized by a low star formation efficiency 
{v ^ OAGyr-^). 

The star formation rate ip{t) is directly proportional the gas fraction G{t) at the time t, 
according to the Schmidt law: 

V'(i) = iyG{t). (2.1) 

Supernovae (SNe) la and II are responsible for the onset of a galactic wind, when the 
thermal energy of the ISM exceeds its binding energy, which is related to the presence of 
a dark matter halo in which the galaxy is embedded (for more details, see Bradamante 
et al. 1998, Lanfranchi & Matteucci 2003). The binding energy of the gas is influenced by 
assumptions concerning the presence and distribution of dark matter (Matteucci 1992). 
A diffuse (i?e/-Rd=0.1, where Re is the effective radius of the galaxy and Rd is the radius 
of the dark matter core) and massive {Mdark/MLum = 10) dark halo has been assumed 
for each galaxy. The time at which the wind develops depends on the assumed star 
formation efficiency. In general, the higher the star formation (SF) efficiency, the earlier 
the wind develops. The models used throughout this paper is characterized by very low 
star formation efficiencies and by a young age, much lower than the times of occurrence 
of the galactic winds. For these reasons, the galactic winds have no effect on the main 
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Figure 1. Observed abundance ratios versus metallicity for the host galaxy of GRB 050730 as 
reported by P07 (solid squares with error bars). The solid line, dashed line and dotted line are 
predictions computed by means of a chemical evolution model for an elliptical galaxy, a spiral 
galaxy and a dwarf irregular galaxy, respectively, and without including dust depletion (see text 
for further details). 
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results obtained in this paper. 

Chemical enrichment from various types of stars is properly taken into account. The 
stellar yields are mainly from Woosley & Weaver (1995) for massive stars, from Meynet 
& Maeder (2002) for low and intermediate mass stars (LIMS) and from Iwamoto et al. 
(1999) for Type la SNe. 

We assume a Salpeter initial mass function (IMF). We assume a cosmological model 
characterized by ^Im = 0.3, J^a = 0.7 and a Hubble constant Hq — 70kms~^ Mpc~^ . 

2.1. Dust production and destruction 

Two elements studied in these paper are refractory: Ni and Fe. To model their gas phase 
abundances, we need to take into account dust production and desctruction. 
The model for dust evolution used in this paper is described in detail in CPM08. 
Here, we summarize its main features. For the refractory chemical element labeled i, 
a fraction Sf^, and is incorporated in dust grains by low and intermediate 
mass stars, type la SNe, and type II SNe, respectively. These quantities are the dust 
condensation efficiencies of the element i in various stellar objects. Here we assume 
Sf^ — 5l°- — = Si — 0.1. This choice is motivated by recent mid-infrared ob- 
servations of one Supernova (Sugerman et al. 2003), which provided an upper limit of 
^ 0.12. The value assumed here is further supported by theoretical studies of the 
local dust cycle (Edmunds 2001) 

Dust grains are usually destroyed by the propagation of SN shock waves in the warm/ionised 
interstellar medium (Jones et al. 1994). If Gdust,i is the fraction of the element i locked 
into dust and G is the gas fraction, the destruction rate is calculated as Gdust,i/Tdestr, 
where Tdestr is the dust destruction timescale, calculated as: 

G 

Tdestr = [eMsMR)-^- • (2.2) 

HSN 

(McKee 1989). Msnr — 13OOM0 is the mass of the interstellar gas swept up by the SN 
remnant. Rsn is the total SNe rate, i.e. the sum of the rates of Type la and Type II 
SNe. 

Unless otherwise specified, we assume that no dust accretion, occurring mainly in cold 
molecular clouds, is taking place in the GRB host galaxy. Our choice is motivated by the 
fact that very little molecular H is observed in local dwarfs, with molecular-to- atomic gas 
fractions of ~ 10% or lower (Clayton et al. 1996). Our assumption is further supported 
by the lack of H2 absorption lines observed in the spectra of GRB afterglows (Whalen et 
al. 2008). In the presence of intense SF, molecular clouds are likely to rapidly dissolve, 
allowing very little dust accretion to occurr. 



2.2. Results and discussion 

The abundance ratios between two elements formed on different timescales can be used 
as " cosmic clocks" and provide us with information on the roles of LIMS and SNe in the 
enrichment process (Matteucci 2001). In particular, the study of abundance ratios such 
as [a/Fe] and [N/a][|Jis quite helpful, since the a-elements (O, S, Mg) are produced on 
short timescales by Type II SNe, whereas the Fe-peak elements and nitrogen are produced 
on Ion timescales by Type la SNe and low and intermediate-mass stars, respectively. As 

f All the abundances between two different elements X and Y are expressed as 
[X/y] — log{X/Y) — log{X/Y)Q, where (X/Y) and {X/Y)q are the ratios between the mass 
fractions of X and Y in the ISM and in the sun, respectively. We use the set of solar abundances 
as determined by Grevesse et al. (2007). 
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Figure 2. Observed abundance ratios versus metallicity for the host galaxy of GRB 050730 
as derived by P07 (sohd squares with error bars) . The triangles stand for upper or lower limits. 
The thick lines represent models with increasing SF efhciency, with the lowest curves having the 
lowest SF efhciency values. The thin dashed lines do not include dust depletion. The red box 
(N/S vs S/H) is used for abundance ratios between non-refractory elements. 

shown by Dessauges-Zavadsky at al. (2004; 2007, hereinafter D07), the simultaneous 
study of the abundance ratios between different elements as functions of a metallicity 
tracer (such as [S/H] or [Fe/H] ) may be used to constrain the star formation history of 
a given system, whereas the study of the abundance ratios versus redshift can be used 
to constrain the age of the system. 

In Fig. 1, we show the predicted evolution of [N/Fe], [Ni/Fc] and [S/Fe] as a function of 
various metallicity tracers, such as [Fe/H] or [S/H], for various star formation histories 
describing galaxies of different morphological types: an elliptical galaxy (solid lines), a 
spiral galaxy (more precisely, a model for the Solar Neigbourhood, dashed lines, see 
CPM08 for a detailed model description) and a dwarf irregular galaxy (dotted lines), 
compared with the values observed for GRB 050730 DLA (solid squares with error bars). 
Fig 1 shows clearly that the most likely progenitor for the host galaxy of GRB 050730 is 
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Figure 3. Observed abundance ratios versus redsiiift for tiie liost galaxy of GRB 050730 as 
derived by P07 (solid squares with error bars). The triangles stand for upper or lower limits. 
The solid (dashed) lines represent a redshift of formation Zf = 6 {zf — 5). The thin dashed 
lines do not include dust depletion. The red box (N/S vs redshift) is used for abundance ratios 
between non-refractory elements. 

a dwarf irregular galaxy, and how the elliptical and spiral model appear indadequate to 
describe the abundance pattern observed by POT in the GRB 050730 host. 
In the following, we will use the dwarf irregular galaxy model to constrain the parameters 
of the SF history and the age of the host galaxy GRB 050730. 

The next step is to constrain the main star formation history parameter of the host 
galaxy of GRB 050730, i.e. its star formation efficiency ly. We search for the model which 
reproduces at best the observed abundances. In figure [2l we show the predicted evolution 
of several abundance ratios for various models, characterized by different SF eficiencies. 
The model which reproduces at best the observed abundances is the one characterized 
by a SF efficiency = 0.01Gyr~^ . 

Once we have determined the star formation efficiency of the best model, the following 
step is to constrain the age of the GRB host galaxy. To perform this task, we study the 



124 



Francesco Calura 



I 

>. 
O 

\ 
K 

CO 

o 




-0.5 



Figure 4. Predicted evolution of tlie SSFR vs [Fe/H] for tiie best model for the host galaxy of 
GRB 050730. The open squares are the values determined for 5 QSO DLAs in the sample by 
D07. 



abundance ratios vs redshift (Fig. [3]) . The abundances are reproduces at best at an ages 
of ~ 0.4 Gyr, corresponding to a redshift of formation Zf ~ 5. At this age, the specific 
star formation rate is ~ 5 Gyr^^, the Dust-to-Gas ratio is ^ 10~^, in agreement with 
the upper limit of 8 • 10^^ as derived observationally by POT. The Dust-to-mctals ratio 
is 0.04, consistent with the upper limit of 0.1 measured by P07. 

Finally, in Fig. 21 we compare the predicted SSFR vs [Fe/H] obtained for the best model 
for the host galaxy of GRB 050730, with SSFR derived for a sample of QSO DLAs by 
D07. The evolution of the SSFR of the GRB host is compatible with the values found 
for QSO DLAs, possibly indicating similar chemical evolution paths. However, this may 
be the result of a coincidence, since in principle GRB DLAs are expected to represent 
denser and more metal-enriched galactic regions than QSO DLAs. By extending the 
method presented here to other systems, it will be possible to shed light on possible 
analogies and differences between QSO and GRB DLAs. 
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